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We present the first direct search for single top quark production using reconstructed tau leptons 
in the final state. The search is based on 4.8 fb _1 of integrated luminosity collected in pp collisions 
at ^/s—1.96 TeV with the DO detector at the Fermilab Tevatron Collider. We select events with a 
final state including an isolated tau lepton, missing transverse energy, two or three jets, one or two 
of them being identified as b quark jet. We use a multivariate technique to discriminate signal from 
background. The number of events observed in data in this final state is consistent with the signal 
plus background expectation. We set in the tau+jets channel an upper limit on the single top quark 
cross section of 7.3 pb at the 95% C.L. This measurement allows a gain of 4% in expected sensitivity 
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channels already published by the DO collaboration with 2.3 fb _1 of data. We measure a combined 
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INTRODUCTION singly by the electroweak interaction. Single top quark 

production can be used to directly measure the CKM 
At the Fermilab Tevatron Collider, top quarks can matrix clement |V th | Q , to determine the top quark par- 
be produced either in pairs by the strong interaction or 



4 



tial decay width and lifetime [2j, to study top quark 
polarization and to probe physics beyond the standard 
model (SM) Q. The production of a single top quark 
is accompanied by a 6 quark in the s-channel mode or 
by both a b quark and a light quark in the i-channel 
mode, as illustrated in Fig. [TJ Besides the s- and t- 



(a) (b) 




dominant background source is multijct events, which 
are poorly modeled by Monte Carlo (MC) simulations, we 
build most of our background model from an independent 
sample of multijct data. We then model several smaller 
background sources using MC and combine them with 
the multijet sample to complete the background model. 
We then train a multivariate discriminant to separate the 
simulated single top signal from the background model. 
Finally, we extract the single top cross section and com- 
bine the result with the existing electron and muon chan- 
nel measurements. 



OBJECT IDENTIFICATION AND EVENT 
SELECTION 



Figure 1: Diagrams for single top quark production in (a) s- 
channel and (b) i-channel production showing the top quark 
decays of interest. 

channels, single top quarks can also be produced in an 
associated tW process via bg — > tW. At the Tevat- 
ron, this channel has a negligible cross section compared 
to s- and t-channel production The s-channel pro- 
cess is referred to as "tb" production, where tb includes 
both tb and tb states. The i-channel process is abbre- 
viated as "tqb" , where this includes tqb, tqb, iqb, and 
iqb states. Considering the SM decay modes of the top 
quark and W boson, single top production and decay 
results in four channels: electron+jets, muon+jets, tau 
lepton+jets (tau+jets), and all-jets. Evidence yHZl and 
observation 0, Q of single top quark production in the 
electron+jets and muon+jets channels and the first dir- 
ect measurement of | V*5 1 [a] have been published recently. 
However, the tau+jets channel has not been measured 
so far due to the overwhelming jet background at the 
Tcvatron, although signatures involving tau leptons have 
been explored by DO, for example, in the measurement 
of the Z — > tt cross section II Ofl and, more recently, in 
the context of Higgs searches The measurement of 
the single top quark cross section in the electron and 
muon channels is still limited by statistical uncertain- 
ties. In this analysis, adding the tau+jets channel in- 
creases the signal acceptance by 32% compared to the 
DO observation In addition, the tau+jets channel is 
a statistically independent channel with different domin- 
ant backgrounds and different systematic contributions 
compared to the electron and muon channels. As such, 
the tau channel provides an independent measurement of 
the single top production cross section. In addition, the 
approach developed in the analysis could be extended to 
other studies, such as Higgs searches in tau+jets chan- 
nels, where the cross section to be measured is low and 
multijet background is dominant. 

In this Letter, we report the first direct search for 
single top quarks in the tau+jets channel. Since the 



The upgraded DO detector is described in detail in 
Ref. [l2|. A right-handed coordinate system is used in 
the analysis. In the system, the z-axis is along the proton 
direction, is the azimuthal angle, r\ is the pseudorapid- 
ity, — In [tan 0/2], where 9 is the polar angle, and the true 
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rapidity is defined as, 1/2 In [(£7 +p z c)/(E — p z c) 
This analysis is based on a sample of DO Run II data col- 
lected between August 2002 and April 2009. Run Ha and 
Run lib data are defined as two sub-datasets correspond- 
ing to integrated luminosities of 1.0 fb _1 and 3.8 fb _1 re- 
spectively. A new inner layer of silicon microstrip track- 
ing detectors was added to the detector between Run Ha 
and Run lib. The additional tracking detectors and the 
increased instantaneous luminosity in Run lib change the 
&-jet identification performance. 

The sample considered contains events which have 
passed one of a list of specialized trigger conditions. The 
most important ones either set a threshold on the total 
scalar sum of transverse momenta of the jets in the event, 
require a minimum transverse momentum of all jets, or 
select events based on the acoplanarity of the two leading 
jets sorted in transverse energy. The trigger efficiency in 
this analysis is «45%. 

A hadronically decaying tau lepton appears as a nar- 
row jet in the DO detector. A tau candidate is a 
calorimeter cluster reconstructed from all the towers 
with energy above a threshold within a cone 1Z = 
\J (A</>) 2 + (A77) 2 <0.5 (where <f> is the azimuthal angle 
and 77 is the pseudorapidity) around a seed tower. The 
tau candidate must have at least one track associated 
with the cluster, and possibly an additional energetic 
subclustcr of cells in the electromagnetic (EM) section 
of the calorimeter [l(|. Hadronic tau candidates are 
separated in three types according to the tracking and 
EM calorimeter information: (1) single track with no 
EM subclusters, (2) single track with EM subclusters, 
and (3) two or three associated tracks. The classifica- 
tion is motivated by the decay modes: (1) r ± — > n ± iy 



(2) 



p ± v and (3) r 



7r ± 7r ± 7r T (7r°)i/. We re- 



quire the tau transverse momentum, pTp, to be larger 
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than 10, 5, 10 GeV for Type 1, 2 and 3 tau lcptons. 
We also require the transverse momentum of the asso- 
ciated track, p^ k , to be larger than 7 GeV (5 GeV) for 
Type 1 (2) tau leptons. For Type 3, the transverse mo- 
mentum of at least one track, p^ k , has to be larger than 
5 GeV and the sum of the associated track transverse mo- 
menta, X)trkPr k ' nas t° be greater than 7 GeV. Hadron- 
ically decaying tau leptons are distinguished from other 
types of jets using variables such as isolation (defined as 
- £^ oro ) I E!j? IC where E T T and £^ oro are the trans- 
verse energy in a cone with radius 1Z = 0.5 and a smaller 
cone with radius 1Z = 0.3 defined about the same axis), 
shower width, and shower profiles (a ratio of the Et's 
of the two most energetic calorimeter towers with size 
Acj) x A?] = 0.1 x 0.1 over Et of the tau candidate). As 
there is no single tau identification variable which can 
provide the required background rejection, a multivari- 
ate technique is used to combine these features into a 
single discriminant. Tau identification is then performed 
by applying kinematic selections as well as a requirement 
on the multivariate discriminant output. 

All other analyses at DO that use hadronic tau de- 
cays except the one reported in this Letter use tau iden- 
tification relying on a neural network (NN) trained on 
Z — > tt decays and background samples suitable for that 
signal [X0| . In contrast, the multivariate technique used 
in this analysis for tau identification relies on boosted 
decision trees (BDT). The BDT technique has been used 
in previous DO single top quark analyses [5, 8j and is de- 
scribed in Ref. @ . In brief, a decision tree is an algorithm 
which combines selection requirements on a large num- 
ber of variables with varying discriminating power into a 
single, more powerful, multivariate discriminant [l3lll4j|. 
It can be "boosted" by building the multivariate discrim- 
inant through a weighted average score from many de- 
cision trees instead of a single tree A total of 25 
well-modeled kinematic variables for each tau type serve 
as the inputs to BDTs. Table U shows the 10 most dis- 
criminative variables with their normalized importance 
values for tau Types 1, 2 and 3. The importance is de- 
rived by an algorithm in which variable usage frequency, 
separation gains and numbers of events in the splitting 
nodes are considered [l6| • A set of trees is created based 
on a simulated tau sample from single top quark MC 
events, and realistic background strongly dominated by 
fake tau leptons. This fake tau background is extrac- 
ted from data by requiring events to pass tau jet triggers 
and applying the kinematic selections given above. Both 
the signal and background have different kinematics from 
the standard NN training samples. By changing the tech- 
nique from neural networks to boosted decision trees we 
gain «3%, «8% and «2% (for Types 1, 2 and 3) signal 
efficiency for the same background rejection rate (98%). 
By changing both the technique and the signal and back- 
ground samples to match the busy single top environment 
with extra jets we gain «8%, ~20% and «8% (for Types 



1, 2 and 3) signal efficiency yielding «76%, ~69% and 
~59% for the same rejection. We require exactly one tau 
lepton per event. 



Jets are reconstructed by an iterative cone algorithm 
with radius 1Z = 0.5 in rapidity-azimuth space [17] . The 
highest-pT jet must have px > 25 GeV and the second 
highest-pr jet px > 20 GeV while any additional jet must 
have px > 15 GeV. The highest-pr jet must have pseu- 
dorapidity |?/| < 2.5 and all other jets |7y| < 3.4. The 
jets and the tau lepton must be isolated by requiring for 
their spatial separation in pseudorapidity-azimuth space 
be larger than 0.5. In order to identify b jets, a neural 
network is trained on the outputs of three 6-jet identi- 
fication algorithms: secondary vertex, jet lifetime prob- 
ability, and counting signed impact parameter [l8j]. All 
three of these algorithms discriminate b jets from light 
quark jets by exploiting the signatures of the relatively 
long lifetime of b hadrons. If the neural network output 
of a jet is larger than 0.775, the jet is tagged as a b jet. 
This operating point corresponds in our selected sample 
to a 6-tagging efficiency of 40% and a light-quark tagging 
rate of 0.4%. We select events with two or three jets, in- 
cluding at least one b jet, in order to enhance the signal- 
to-background ratio. We also require 20 <-$^< 200 GeV 
where $t is the missing transverse energy which is equal 
to the negative of the vectorial sum of the transverse 
energy deposited in the calorimeter by all particles. A 
tau energy scale correction has been applied and $ T has 
been corrected for the presence of the tau leptons. We do 
not exclude electrons that satisfy the tau identification 
requirement since these electron events provide >50% of 
our signal acceptance. However, we veto events with one 
isolated electron or one isolated muon to make sure the 
tau+jets sample has no overlap with the electron and 
muon samples in order to be able to combine the meas- 
urements. The data have been split by tau (Types 1 and 
2 combined and Type 3), jet multiplicity (two jets and 
three jets), number of b jets (one b jet and two b jets) 
and running period, for a total of 16 analysis channels. 

We select 3845 6-tagged tau+jets candidate events, 
among which we expect 72 single top quark events. 
Table im shows the event yields for all channels combined. 
About 85% of single top quark events in this sample come 
from tau Types 1 and 2 and 86% are events with only 
one b jet. The acceptance times efficiency is 3.0% when 
considering only hadronic tau leptons. 



SIGNAL AND BACKGROUND MODELING 

Single top quark events are simulated by the next- 
to- leading order (NLO) event generator SINGLETOP fljj] . 
which is based on COMPHEP 20l.[21|. 
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Table I: The 10 most discriminative variables with their normalized importance values in the training of the tau identification 
BDT. The variables listed are explained in Appendix I. 



Rank 


Tau Type 1 


Importance 


Tau Type 2 


Importance 


Tau Type 3 


Importance 


1 


Widths (t) 


2.1X10" 1 


Widths (t) 


2.8 xlO -1 


Isolation 


7-lxlO- 1 


2 


Isolation 


2.0X10" 1 


Profile 


2.8X10- 1 


Rati0 T; trks 


3.9x10-2 


3 


Ratio T: trks 


9.1xl0" 2 


ei2 


l.SxlO" 1 


ei2 


3.0x10-2 


4 




6.9X10" 2 


Profile' 


5.5xl0" 2 


TTialltrk— trkl — 

H/rp 


tTk2 /E^ 2.5x10-2 


5 


Profile' 


5.1xl0- 2 


5a 


2.3xl0- 2 


_trkl 
Z DCA 


1.9x10-2 


6 


£,EMlaycr3 


3.5xl0~ 2 


Profilei a yer3 


2.3xl0- 2 


Widths (t) 


1.6x10-2 


7 


Isolation' 


3.2X10" 2 


Ratio T) t x ks 


2.0x10^2 


77'trkl / rpr 

Ulrp 1 Hjrp 


1.5x10-2 


8 


Profile 


3.2xl0- 2 


77iEMcl2 / TPT 

Jl/rp j lljji 


1.7x10-2 


Profilei aycr3 


1.3x10-2 


9 


RatiOEM12,r 


3.0xl0" 2 




1.7x10-2 


Width^(r) 


1.1x10-2 


10 


Width'^(r) 


2.8xl0" 2 


Isolation 


1.5x10-2 




1.1x10-2 



Table II: Expected and observed events in 4.8 ft)- 1 of integ- 
rated luminosity shown in tau Types 1 and 2, Type 3 channels 
and all analysis combined. The uncertainties include both 
statistical and systematic components. 



Source Types 1 and 2 Type 3 Sum 

tb+tqb 61 ± 11 11 ± 2 72 ± 12 

VF+jcts 573 ± 68 107 ± 12 680 ± 104 

Z+jets 43 ± 8 17 ± 3 60 ± 10 

Dibosons 30 ± 5 7 ± 1 37 ± 6 

U 170 ± 35 60 ± 12 230 ± 44 

Multijcts 1444 ± 38 1182 ± 21 2626 ± 98 

Total prediction 2321 ± 94 1384 ± 28 3705 ± 153 

Data 2372 1473 3845 



Since tau leptons arc observed as narrow jets of 
particles in the calorimeter, the main background to 
single top quark events in the tau+jets channel is mul- 
tijet production. This is unlike the other leptonic single 
top channels in which H^+jets events are the main back- 
ground [1, IH . We have developed a method to model the 
multijet background directly from data. The principal 
steps in this method can be summarized as: 

1. Derive a tag rate function (TRF) to describe 
the probability to 6-tag any individual jet in the 
sample. 

2. Apply this TRF to the data sample that has no 
^-tagged jets. 

3. Using simulated events for other physics sources, 
subtract them from the sample derived in Step 2 
to get "purc-multijcts" . 

4. Normalize the sample derived in Step 3 to data. 

5. Combine the derived background sample, pure- 
multijets, with simulations of other background 



sources: tt, T-F+jets, Z+jets, dibosons. 

In Step 1, we take the ratio of the number of 6-tagged 
jets in our data sample to the total number of jets to 
define a tag rate: the average probability that a jet is 
identified as a b jet. We measure the tag rate as a function 
of jet pt and r\ and jet multiplicity. 

In Step 2, we apply these TRFs to those events that 
have no 6-taggcd jets. This TRFcd sample is kinematic- 
ally similar to our analysis sample, but there is no overlap 
since we require at least one 6-tagged jet in our analysis 
sample. 

In Step 3, we remove physics background sources such 
as ti, W+jets, Z+jets and dibosons. In this procedure, 
we subtract from the zero-tagged TRFed multijet sample 
the contaminations of ti, VF+jets, Z+jets and dibosons. 
Other background sources are modeled through simu- 
lations. These simulations, except the tau decay, have 
been described in Q . The program TAUOLA [22| (version 
2.5) was used to model the decays of tau leptons includ- 
ing polarization effects. We normalize the TU+jcts back- 
ground to match data by the scale factors that are derived 
from the study in the clcctron+jcts and muon+jcts chan- 
nels [H . We apply TRFs to the zero-tagged MC samples 
to estimate the contamination mentioned above. A sim- 
ilar procedure is used to ensure that any small single top 
signal contamination in the background data sample is 
also subtracted. 

In Step 4, the multijet events after contamination re- 
moval are normalized to data in a multijets-enriched re- 
gion, as defined by the background-dominated region of 
the multivariate discriminant described below. 

In addition to multijet events modeled by the proced- 
ure described above, our background model includes tt, 
W+jets, Z+jets and dibosons modeled directly from sim- 
ulation. In Step 5, we combine these simulated samples 
with the data-derived multijet sample. 

At the end of the background modeling procedure, 
we investigate approximately 150 topological variables to 
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confirm that data and the background model are in good 
agreement since it is expected that the single top quark 
events represent only a small fraction, «2%, of the selec- 
ted data sample. The variables can be categorized in four 
classes: object kinematics, jet reconstruction, top quark 
reconstruction and angular correlations. Figure [2] shows 
four discriminating variables: W boson transverse mass, 
tau transverse momentum, azimuthal angle between the 
second- highest-pr jet and $ T , and cosine of the angle 
between the tau lepton and a jet candidate that is used to 
reconstruct the best top quark mass (defined as closest to 
170 GeV). These variables are shown for the most sens- 
itive channel: Types 1 and 2, two jets, one of them b 
tagged. 

BOOSTED DECISION TREES 

It is expected that single top quark events are only 
a small fraction of the selected data sample. We use 
the BDT technique to separate the signal from the back- 
ground. We also employ the Kolmogorov-Smirnov (KS) 
test [23| to verify the compatibility of discriminating 
variables in shape between data and background model. 
From the sal50 variables studied, 44 to 70 of them are 
selected as input variables to train BDTs depending on 
the individual analysis channels. We select only variables 
which have a KS probability > 0.1. The KS values of the 
selected variables are uniformly distributed above this 
value. 

Separate sets of BDTs are built with these variables for 
each analysis channel. Table inil lists the 15 most discrim- 
inative variables with their normalized importance values 
in the most sensitive channel. Figure [3] shows the BDT 
output with all channels combined, in the region between 
0.5 to 1.0, i.e. where the single top quark signal events are 
expected. Data and the background model are in good 
agreement in the region. The background-dominated re- 
gion from 0.0-0.2 is used to define the multijets-enhanced 
region used in Step 4 of the multijet background model- 
ing procedure. 

SYSTEMATIC UNCERTAINTIES 

We consider systematic uncertainties from correction 
factors applied when modeling the signal and back- 
ground [6|. "Normalization" uncertainty components 
from the correction factors affect the signal efficiency 
and the normalization of the background samples, while 
"shape" uncertainties change the shapes of the distribu- 
tions for the background and the expected signal. The 
largest uncertainties arise from VF+jets normalization to 
data, tau identification efficiency, tag rate functions, and 
jet-flavor correction in IF+jets and Z+jets events. Other 
uncertainties include multijets normalization, integrated 



Table III: The 15 most discriminative BDT training variables 
with their normalized importance values in the most sensitive 
channel. A</>(objl, obj2) is the azimuthal angle between objl 
and obj2. cosa(objl, obj2) is cosine of the angle between 
objl and obj2. "jetl" and "jet2" are the highest-pr jet and 
the second- highest-pr jet, respectively. "jetl+jet2" is a sys- 
tem consisting of "jetl" and "jet2". The subscript, "Top- 
Frame", indicates that the reference frame is the rest frame 
of a top quark which is reconstructed using a b-tagged jet, 
while the subscript "tag" ( "untag" ) refers to the jet passing 
(failing) the b-jet identification algorithm. Xltrks in C vtPT k is 
the transverse momentum of the vectorial sum of all tracks 
with a cut on the distance of closest approach (DCA) to the 
primary vertex. 



Rank Variable Importance 

1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 



luminosity, MC statistics, jet energy scale, jet identific- 
ation, jet energy resolution, initial- and final-state radi- 
ation, jet fragmentation, theoretical cross sections, the 
rewcighting of the jet angular distributions in IF+jets 
events, signal contamination removal, non-multijets con- 
tamination, branching fractions, instantaneous luminos- 
ity reweighting, parton distribution functions, primary 
vertex selection, and tau energy scale. The total uncer- 
tainty on the background model is 4.2%-19% depending 
on the analysis channel. Table llVl summarizes all sources 
of uncertainties considered. 

Some of the uncertainties are common with the study 
in the electron and muon channels and have been presen- 
ted in Ref. @ . Below arc the uncertainties specific to this 
analysis: 

(i) W +jets normalization to data (normalization) 
The uncertainty is on the scale factors applied to 
normalize VF+jets to match data. Since we use 
the scale factors derived from the electron+jets and 
muon+jcts study we consider the difference 
between these two channels as the uncertainty in 



W boson transverse mass 


5.0x10" 


-1 


A<Kjet2, $ T ) 


3.0x10" 


-1 


Pt(t) 


3.3x10" 


-2 


cos a(best jet, r) 


2.8x10" 


-2 


p T (jetl+jet2) 


1.8x10" 


-2 




1.5x10" 


-2 


COSa(tag, r)To P Framo 


1.1x10" 


-2 


A0(jetl, $ T ) 


1.1x10" 


-2 


Z^trka in evt PT 


1.0x10" 


-2 


Best top quark mass 


1.0x10" 


-2 


PT(best jet) 


7.6x10" 


-3 


Q(t) x r;(untag) 


5.6x10" 


-3 


z position of primary vertex 


4.9x10" 


-3 


K{t, jet2) 


4.6x10" 


-3 


TZ(t, jetl) 


3.9x10" 


-3 
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(c) 
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cosa(BestJet,x) 



(d) 



Figure 2: Comparison between data and background distributions in the most sensitive channel: two jets, one b tag, tau Types 
1 and 2 combined, (a) W boson transverse mass (b) tau transverse momentum (c) azimuthal angle between the second-highest- 
Pt jet and $ T , and (d) cosine of the angle between the tau and the best jet candidate that is used to reconstruct the best top 
quark mass (defined as closest to 170 GeV). In (b), the double-peak structure is caused by different pr threshold for tau types 
1 and 2. 



the tau+jets channel. 

(ii) Multijets normalization (normalization) 

The statistical uncertainty of the multijet sample in 
the BDT region [0.0, 0.2] is used. 

(hi) Tag rate functions (shape and normalization) 

This uncertainty consists of two components: those 
on the multijet background sample and those on 
the MC samples related to 6-tag modeling. The 
former is evaluated by raising and lowering the tag 
rate by one standard deviation of its experimental 
determination. Uncertainties considered in the lat- 
ter are from several sources: statistics of the simu- 
lated events; the assumed heavy flavor fractions in 
the simulated multijet sample used for the mistag 
rate determination; and the choice of paramcteriz- 



ations @. 

(iv) Tau identification efficiency (normalization) 

This uncertainty is estimated by the difference in 
tau identification efficiency between data and MC 
as derived in a tau-enriched data sample. 

(v) Signal contamination removal (shape and normal- 
ization) 

In Step 3 of the background modeling, we reweight 
single top quark events to remove any small sig- 
nal contamination. The uncertainty is evaluated 
by raising and lowering the weighting function by 
one standard deviation. 

(vi) Non-multijets contamination removal (shape) 

In Step 3 of the background modeling, we subtract 
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tb+tqb boosted decision tree output 

Figure 3: Distribution of BDT output with all channels com- 
bined in the signal region (BDT>0.5). The single top quark 
signal (tb + tqb — »tau+jets) is normalized to the measured 
cross section. 

the non-multijets contamination from the zero- 
tagged TRFed multijet sample by weighting events. 
The uncertainty is evaluated by raising and lowering 
the weighting function by one standard deviation. 

(vii) Tau energy scale (normalization) 

The energy of hadronic tau candidates with low 
energy is corrected using the energy in the calor- 
imeter and the momentum of the tracks associated 
to the tau leptons using parameterized single pion 
response functions. The uncertainty on the scale is 
estimated by varying these parameterizations. 



RESULTS 

The number of events observed in data and the shape 
of the BDT discriminant are consistent with the sum of 
the signal and background predictions. To estimate the 
statistical significance of the signal observation we use 
the same Bayesian approach as in Refs. H Q. This in- 
volves forming a binned likelihood as a product over all 
bins and channels. When measuring a cross section, its 
central value is defined by the position of the peak in the 
posterior density, and the 68% interval about the peak is 
taken as the uncertainty. The posterior density is integ- 
rated from until 95% of the posterior area is contained 
and the upper limit is set at this point. Systematic un- 
certainties, including all correlations, are reflected in this 
posterior interval. Assuming a single top quark cross sec- 
tion of 3.46 pb for a top quark mass of 170 GeV [3], we 
estimate the expected sensitivity to the standard model 
signal by calculating the ratio of the position of the peak 
of the expected posterior density to its lower half width. 



Table IV: A summary of the relative systematic uncertain- 
ties for each of the correction factors or normalizations. The 
uncertainty shown is the relative error on the correction or 
the efficiency, before it has been applied to the MC or data 
samples. We do not show relative systematic uncertainties of 
the components for shape since they depend on distribution 
binning. 



Relative Systematic Uncertainties 


Components for Normalization 




b-jet fragmentation 


2.0% 


Branching fractions 


1.5% 


Diboson cross sections 


5.8% 


Instantaneous luminosity reweighting 


1.0% 


Integrated luminosity 


6.1% 


Initial- and final-state radiation 


(0.6-8.0)% 


Jet energy resolution 


4.0% 


Jet energy scale 


(4.0-14.0)% 


Jet fragmentation 


5.0% 


Jet identification 


1.0% 


MC statistics 


(0.5-16.0)% 


Parton distribution functions 


3.0% 


(signal acceptances only) 




Primary vertex selection 


1.4% 


Multijets normalization 


(3.0-14.0)% 


Tau energy scale 


(1.0-1.5)% 


Tau identification efficiency 


11.0% 


Triggers 


5.5% 


ti cross section 


12.7% 


IF+jets heavy-flavor fraction 


13.7% 


VK+jets normalization to data 


(7.0-15.0)% 


Z+jets cross section 


3.6% 


Z+jets heavy-flavor fraction 


13.7% 


Components for Shape 




Alpgen reweighting on W+jets sample 




Non-multijets contamination removal 




Components for Shape and Normaliza- 




tion 




Signal contamination removal 




Tag rate functions 





This yields a ratio of 1.8, i.e. a sensitivity corresponding 
to approximately 1.8 standard deviation. 

In order to test the linearity of our procedure with 
respect to the single top quark cross section, we gen- 
erate several ensembles of pseudodatasets by randomly 
sampling from background model events. We specify five 
input signal cross sections: 2.0 pb, 3.46 pb, 6.0 pb, 8.0 pb 
and 10.0 pb and generate ensembles at each value. Each 
ensemble contains w2000 pseudodatasets with all system- 
atic uncertainties considered. We then measure the cross 



10 



-0.22F" 



DO Run II 4.8 fb 

""pb 

95% C.L. Upper Limit=7.3 pb 

3.7 ± 2.2 pb 
95% C.L. Upper Limit=8.5 pb 




Signal Cross Section [pb] 

Figure 4: Expected SM and measured Bayesian posterior 
probability densities for the tb+tqb cross section. The shaded 
regions illustrate ± one standard deviation from the peak loc- 
ations. 



section in each of the 2000 pseudodatasets at each input 
value and assess linearity. A linear fit to the measured 
vs. input cross section gives a slope of 0.99±0.01 and 
intercept of —0.14 ±0.05. Therefore, over the range con- 
sidered, there is no significant evidence of bias in the 
measurement procedure. 

We obtain an observed posterior density that is used 
to define an upper limit on the cross section assuming 
no signal. We can use the same technique to determine 
an observed cross section and its uncertainty. Assuming 
no signal, we extract an upper limit of 7.3 pb at 95% 
C.L. If we perform a cross section measurement, we ob- 
tain 3.4^'g pb. The measured sensitivity, a ratio of the 
position of the peak of the measured posterior density 
to its lower half width, is 1.9. Figure [4] shows the expec- 
ted and measured posterior densities with shaded regions 
corresponding to ± one standard deviation from the peak 
locations. 



COMBINATION WITH OTHER CHANNELS 

As this data sample has no overlap with that used 
in H, it is straightforward to combine the results. 
In the combination, the tau channel and the (elec- 
tron, muon)+jets channels are treated as two independ- 
ent channels using the same Bayesian approach used to 
combine different tau channels above. The ratio of the 
position of the peak of the expected posterior density to 
its lower half width is 4.7, compared to 4.5 in the elec- 
tron+jets and muon+jets channels combined. We gain 
4.4% in expected sensitivity by adding the tau+jets chan- 
nel. The observed posterior density is also calculated and 
yields a combined cross section of: 



a{pp ->tb + X, tqb - 



X) = 3.84™ pb 



Figure [5] shows several recent measurements of single top 
quarkproduction compared to the theoretical SM predic- 
tion [J], 3.46±0.18 pb, calculated for a top quark mass 
of 170 GeV 0. 



DO e+|i combination (2.3 fb"') . « 3.94 tSlSs P b 

DO i+jets (4.8 fb' 1 ) , 3.40 pb 

DO e+u+T combination 
CDF combination (3.2 fb'V 
CDF+D0 combination 

Theoretical SM prediction at top quark mass 170 GeV 



4 



1 2 3 4 5 6 

a(pp-> tb+X, tqb+X) [pb] 



Figure 5: Summary of several recent measurements of single 
top quark production cross section. The theoretical SM pre- 
diction Q at a top quark mass of 170 GeV is included as 
a shaded band. The "DO e+fi combination" result is taken 
from while the "CDF combination" result comes from [s| 
and the "CDF+D0 combination" result from 1261 1. 



tt CROSS CHECK 



As an additional cross check of our background model, 
we have measured the top quark pair production cross 
section in the same data sample, including systematic 
uncertainties and using the same background model and 
the same techniques as we use to measure the single top 
quark cross section. We measure a top quark pair pro- 
duction cross section of 10.01-^6 pb, in good agreement 
with the theoretical expectation, 7.9lj[Q;5g pb, from the 
next-to-next-to-leading-order calculation for a top quark 
mass of 170 GeV 24 \, and a recent DO experimental res- 



ult, 8.18^j]'g7 pb, also for the same top quark mass [25| 



SUMMARY 

In summary, we have presented the first direct study 
for single top quark production in the tau+jets channel 
using 4.8 fb^ 1 of integrated luminosity at the DO experi- 
ment. Due to different dominant backgrounds and differ- 
ent systematic uncertainties from the electron and muon 
channels, the tau+jets channel serves as a channel to in- 
dependently search for single top quarks. To increase 
sensitivity, electron+jets events not entering the meas- 
urement in the electron+jets channel and where the elec- 
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tron satisfies tau identification criteria are also included 
in the tau+jets sample. An upper limit of 7.3 pb at the 
95% C.L. for the cross section is obtained. The expected 
sensitivity of the tau+jets channel alone is 1.8 stand- 
ard deviations. Adding the tau+jets channel increases 
the signal acceptance by 32% compared to the DO ob- 
servation analysis, which was based on electron+jets and 
muon+jets channels. The expected sensitivity of the elec- 
tron+jets, muon+jets and tau+jets combined analysis is 
4.7 standard deviations, to be compared to 4.5 standard 
deviations in electron+jets and muon+jets alone. The 
measured cross section in all three combined channels is 
found to be 3.84±g;|| pb. This is the most precise meas- 
urement to date of the single top quark production cross 
section. 
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Appendix I: Tau Identification Variable Definitions 



Isolation 



x . rk , where ^ Pt ^ s ^ e sum °f Pt °f 
non-tau-associated tracks within a cone size 0.5 and 
p^ rk is the sum over all tau-associated tracks' px- 

Isolation': if |?7 d et| ^ 1-0, where 7yd e t is tau's detector 
pseudorapidity, which is defined with respect to the 
center of the detector, Isolation' = Isolation. If 
\Vdet\ > 1-0, Isolation' = Isolation/ (1.5 x ?7dot| — 
0.5). 

Profile: e ti+e T 2 wnere jr; T1 an d £ T2 are the trans- 
verse energies of the two highest-pT calorimeter 
towers in a tau object. 

Profile': if |r? det | <1.5, Profilc'=Profilc. If |r? det | >l-5, 
Profile' = Profile x (0.67 + 0.22 x \r/ dct \). 

Profilei aycr 3: a ratio of Et of the highest pt EM sub- 
cluster over Et deposited in the 3 rd layer of the 
EM calorimeter within a cone R < 0.5. 



RatiOEM12,r = 



where E KM1 and E EM2 are 



energies deposited in the 1 st and 2 nd layers of the 
EM calorimeter. 



Ratio r trks^ 



Width J?! 0(r): tau shower width, the root sum of squares 
of the -fir-weighted r\-<\> distance of all calori- 
meter towers with respect to the tau axis, i.e., 



Y^i=\ i^Vi + ) ; §p where i is the index of 
calorimeter towers and Et — Xa-^Ti- 

Width^(r): Width^(r)/(1.0 + 0.29x |ry det |). 
2 dca : 2 position of the highest-pr track at DC A. 



5a: \J (Acf>) 2 + (A?y) 2 , where Acf> and Arj arc differences 
between the vector sums of tau tracks and of all 
EM subclusters O. 



e 12 : yE T tA ■ E EMcls /E T , where E T tA is the sum over 

all tau-associated tracks' Et, E™ c1b is Et of the 
sum over EM subclusters. For a system of tau- 
associated tracks and EM subclusters, the observed 
tau transverse mass is ei2 x E T x da in the small 
angle approximation [loj |. 

r-ialltrk — trkl — trk2 , r , i j_ 

hi T : transverse energy of tracks except 

the first two highest-px tracks. 

^EMci2. transverse energy of the second-highest-py 



EM subcluster. 



E, 



EMlayer3. 



transverse energy deposited in the 3 layer 



of the EM calorimeter within a cone 1Z < 0.5. 



E^ kl : transverse energy of the highest-pr track. 
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